disp("ProblemC4:")
for n = 1 : 6
    disp("n=");
    disp(n);
    H = hilb(n+1);
    disp("cond(H)=")
    disp(cond(H));
    alphaP4 = realHilbert(n);
    disp("alphaP4=");
    disp(alphaP4);
end